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Abstract. - The different kinds of behaviour of three-body systems in the weak binding 
limit are classified with specific attention to the transition from a true three-body system to 
an effective two-body system. For weakly bound Borromean systems approaching the limit of 
binding we show that the size-binding energy relation is an almost universal function of the 
three s-wave scattering lengths measured in units of a hyperradial scaling parameter defined as 
a mass weighted average of two-body equivalent square well radii. We explain why three-body 
halos follow this curve and why systems appearing above reveal two-body substructures. Three- 
body quantum halos 2-3 times larger than the limit set by zero hypermoment are possible. 



Introduction. - Attempts for a general classification of halo states were started early in 
the development of the field, see e.g. GJ], and have recently led to the suggested definition of 
halo states as having more than 50% probability of being in a cluster configuration where 
more than 50% of this probability should be in a classically forbidden region. This definition 
is straightforward to apply for a two-body system, where one basically has to find the outer 
classical turning point in the radial motion. The three-body systems are more challenging [Q 
and it is the purpose of this paper to discuss their possible modes of behaviour. 

One obvious point that needs clarification is how to generalize the outer classical turning 
point, that also can be used j^j to scale different halo systems so that e.g. nuclear and 
molecular halos can be compared in dimensionless units. For the interesting special case of 
Efimov states a universal scaling property predicting one Efimov state from the previous one 
has been developed, see Another point that has only been discussed briefly in the 

literature so far is how the transition from a three-body to a two-body state takes place as the 
binding potentials are changed || . Connected to this is the classification of possible three- 
body configurations into Borromean P,^C|,pT|, tango |l2| or other bound states. To clarify the 
principles we shall mainly consider systems where all particles are in relative s-waves, which 
dominate at large distances. Results of more realistic calculations will also be given for 11 Li 
and the hypertriton. 
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Three-body systems. - There is naturally more variability in three-body systems than in 
two-body systems. The three two-body subsystems might all play a role in the asymptotic 
region, so even in the weak-binding limit we should expect several types of behaviour to be 
possible, even for the simplest case of only relative s-waves. Systems with zero and one bound 
subsystem are called Borromean and tango systems, respectively. 

In a two-body system the classical turning points are found by equating the total energy 
and the potential energy. In principle we can generalize this to three particles in a specified 
quantum state described by the wave function ^ and find the probability for being in the 
non-classical region as J |\&| 2 dr, where the integration is confined to regions with potential 
energy larger than the total energy. This will be much harder to calculate than for a two- 
body system and could therefore be a rather impractical condition. Furthermore, we are not 
assured that the wavefunction will behave in a simple way in the non-classical region; there 
might be configurations where one pair is in a forbidden region whereas the third particle is 
in an allowed region. We shall therefore here rather explore the possibility of generalizing 
the classical turning point into a three-body scaling radius and discuss two different types of 
"derivation" of it. 

First derivation. - We use hyperspherical coordinates to describe the relative motion of 
three particles with masses m,-, where i = 1,2, 3. The total mass is M = mi + m-i + m^, the 
individual momenta and coordinates are pi and v, and the hyperradius p is defined by 

mp 2 = i ^2 m r m k (r ; - r k ) 2 = ^2 m i ( r i ~ R ) 2 > C 1 ) 

i<k i 

where R is the center of mass coordinate and m is a mass unit chosen for convenience. The 
hyperradius is an average radius coordinate, applicable to all three-body systems and useful 
for all angular momenta and for non-spherical systems. The total mean square radius (r 2 ) is 
then via the particle sizes (r 2 )i given by 

M(r 2 ) =TO(p 2 )+^™ l (r 2 ) J . (2) 

i 

It is natural to choose a three-body scaling radius po so that the arbitrary mass m enters in 
the same way in p and po and all measures of size, that typically rely on their ratio, become 
independent of m. The two-body scaling property relating size and binding energy can then be 
generalized if (p 2 / p 2 ,} is an almost single- valued function of another dimensionlcss quantities 
Bmp^/Ti 2 , where B is the three-body binding energy ||. Such a scaling property is clearly 
an advantage when searching for a general definition of the scaling radius po- The relation 
should apply for systems consisting of particles with widely different masses and ranges of 
interactions. A tempting definition of po is to maintain the complete analogy to p, i.e. 

i<k 

where Rik is interpreted as the equivalent square well radius of the system consisting of particle 
i and k. As argued in || this definition is convenient in descriptions of three-body systems 
intermediate between two and three-body scaling. 

Hyperradial potential. - The choice of hyperspherical coordinates leads to effective radial 
potentials obtained by adiabatic expansion or by averaging in other ways over the remaining set 
of angular coordinates. The classically allowed regions for such one-dimensional potentials are 
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easily defined. However, they could be completely different from the regions where the three 
pairs of particles are located in their classically allowed regions defined by the corresponding 
two-body potentials. In fact, it is entirely possible to have classical motion in the hyperradial 
coordinate while the system is in non-classical regions in real space. Hyperradial turning points 
are therefore useless as definitions of quantum halos and often without any resemblance to 
the length unit po defined above in terms of two-body properties. 

It is instructive to consider the general behaviour of the hyperradial potential. For zero- 
range two-body potentials the only energy available through combination of parameters is 
ft 2 /(2/ip 2 ), where p is a combination of reduced masses. With this large distance behaviour 
any number of solutions is possible, ranging from zero to the infinitely many Efimov states. For 
finite range interactions the ranges or alternatively the scattering lengths provide additional 
length parameters and the hyperradial potentials could approach zero faster than p~ 2 . 

From |p^| we obtain the large distance behaviour of the dominating lowest s-wave adiabatic 
potential as 





where the reduced mass is pik = miink / (m$ + m/s) 
the average scattering length is defined as 



the s-wave scattering length is a^ and 




(5) 



i<k 



The chosen normalization reduces a av to the common an- when all three particles are identical. 
Thus the large distance behavior of V a d in Eq.(U) is p~ 3 which is reached when p is comparable 



to 48a a „/(7rv / 2), see (jj . 



Second derivation. - Scaling can be shown analytically to occur for the special case of a 
K = wavefunction (K being the hypermoment) for square well two-body potentials of depth 
Sik and radius Rik- Here the effective hyperradial potential V can be obtained Q as 




(6) 



which is valid when p is several times larger than any of the square well radii. The square well 
two-body s-wave scattering length is given by aik/Rik = — 1 + tan (Kn-Rik) / (KitRik) , where 
K~ik is the zero energy wave number inside the square well, i.e K 2 K 2 k /(2pik) = Sik- Then 
K 2 k R 2 k = 2SikPikR 2 k/^ 2 i s a specific function of aik/Rik, which approaches the constant 
7r 2 /4 when ay, becomes much larger than Rik- Thus the effective radial potential in Eq.(^|) 
approaches the form 





where the definition of po, 




(8) 



i<k 



reduces to that of Eq.(||) for identical masses and radii. Note that Eq.(||), in contrast to Eq.(^), 
employs a linear (not squared) summation. For systems where the two-body scattering lengths 



4 



EUROPHYSICS LETTERS 



\ B<B 12 
^^0<B 12 

Tango 

0<B ]2 <B 


/ unbound 
/ B<0 

/ B>0 Borromean 


^Efi 

a 12 = a 13 = a 23/' 


nov ^ l/a 12 

0<B 23 =B ]3 <B 

l/a 23 =l/a 13 



Fig. 1 - Sketch of the different regions of stability for a three-body system as functions of the inverse 
s-wave two-body scattering lengths a^. The central point corresponds to all = oo which is 
assumed to be the threshold for binding of the first state in the i — k subsystem. All potentials are 
attractive or vanishing. Particles 1 and 2 are identical. The dashed line represents three identical 
particles. The two and three-body binding energies are B ik and B, respectively. 

are large, most of the wavefunction will reside in the region where Eq.(Q) holds for K — and 
one clearly has scaling. 

A desirable property of the scaling radius po would be that classically forbidden regions 
on average are given roughly by p > po. We can consider this question briefly for the specific 
K = solutions, where the probability Pik for particles i and k being inside their square well 
radius Rik for large p is 




This probability is 1/2 for p = 2(4/ (in)) 1 / 3, / mRik which indicates that systems reside 
mainly in classically forbidden regions if their mean square radii are somewhat larger than p\. 

We do not believe that the K = wave functions are realistic solutions for the weakly 
bound systems since they do not allow for any form of correlations between the particles. For 
very loosely bound systems where a av is much larger than pq the potential will fall off as p~ 2 
between po and a av , i.e. slower than p~ 3 as for a pure K = solution. 

Scaling properties. The scaling radius will be used to look for scaling properties of 
three-body systems in the weak binding limit. We are aiming for as universal properties as 
possible, but the rather different types of structure that occur in three-body systems means 
we first have to look at how they can be classified. In fig.|l| we illustrate the various stability 
regions as function of scattering lengths. At the border between tango and Borromean regions 
one subsystem has a bound state with zero energy. On the thick horizontal line two identical 
subsystems have bound states of zero energy and the infinitely many Efimov states arise. 

The figure only shows the region where 023 = 013, appropriate e.g. when particles 1 and 2 
are identical, but already indicates that two distinct types of transitions can occur as a function 
of ai2, namely moving from the tango region (the 23 and 13 subsystems are unbound) from 
left to right either directly into unbound systems or through the Borromean region. To see 
this in detail we now turn to the numerical results and show in fig. || the region of weak 
binding for a number of both schematic and realistic examples. Through Eq.(0) this scaling 
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Fig. 2 - Scaling plot for three-body halos. The ratio of the halo and effective potential mean square 
radii is plotted versus the scaled separation energy. The definition of scaling radius in Eq.(^) is 
used. The dashed line represents the Efimov states for minimum attraction. Triangles are for masses 
corresponding to 11 Li. Squares and circles are for ^H. The realistic points are indicated by a large 
closed triangle and square, respectively. The stars and crosses refer to a system of three different 
particles with two fixed scattering lengths while the third is varied. Almost indistinguishable curves 
arise for three identical Borromean bosons by varying the common scattering length. The arrows 
indicate transitions from (i) Borromean to tango region (closed circles) , (ii) tango directly to unbound 
three-body system (open circles), (iii) Borromean to either two bound subsystems or the tango region 
(stars, crosses, open and closed triangles and squares) all occurring approximately at the same point. 
The excited states of 11 Li all correspond to a bound n- 9 Li system. 



plot displays the mean square radius of the system in units of p\ versus the dimensionless 
binding energy mBp^/h 2 where B is the three-body binding energy. Zero binding energy 
corresponds to three non-interacting particles at rest. 

We shall in this letter use the scaling radius defined in Eq.(|§|), where Rik is the radius 
of the square well with the same scattering length and effective range as the actual two- 
body potential. Changing to the other definition in Eq.(^) has essentially no effect on the 
hypertriton examples, whereas the different n Li points move by about 30%. The resulting 
figure is practically indistinguishable from the present one, except that the perfect scaling for 
K = is violated slightly. There is no strong preference for any of these two definitions of po. 

The first analysis in Jj| assumed that simple systems at least aymptotically would corre- 
spond to a single value of the hypermomcnt K. We have artificially restricted the wavefunc- 
tions for different systems to contain only K=0 terms and find that these systems, as argued 
above, indeed do scale and lie on a single curve. However, the points corresponding to more 
realistic calculations lie above this K=0 curve and can be grouped and understood as follows. 

In the realistic ^H-structure the neutron and proton are almost in a deuteron configuration 
(close circles) whereas the A-particle is far outside but still bound by about 0.14 MeV. Thus 
this is a tango system and the large size is almost entirely due to the A-deuteron extension, 
i.e. of two-body character. By using a simplified form for the neutron-proton interaction 
with only s-waves the position is sligthly higher than the realistic point. By decreasing 
this neutron-proton attraction the binding decreases moderately while the radius drastically 
increases (open circles). This off hand surprising property can be understood from fig. [I] by 
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moving horizontally from the tango region to the right, increasing l/a\2- For sufficiently 
weak initial binding, i.e. large values of l/ai3, the threshold for A-deuteron binding is reached 
instead of the Borromean region as we decrease the attraction of the only bound subsystem. 
Thus the diverging radius is due to A-deuteron two-body threshold and not related to the 
neutron-proton threshold. 

Decreasing the s-wave attraction between neutron and proton while maintaining all other 
parts of the realistic interaction (filled circles), now leads from tango into the Borromean region 
of fig-0, i.e. the three-body system remains bound even after the neutron-proton potential 
is too weak to form a bound two-body state. The result is a completely different curve in 
fig. much more in agreement with the logarithmic divergence expected from although 
still significantly above the K = curve. We stress that nothing drastically happens at the 
arrow where the deuteron becomes unbound (scaled binding energy about 0.08). 

The realistic point for the Borromean nucleus n Li is at a scaled binding energy of 0.13 
in fig. |^. By decreasing the s-wave attraction in the neutron- 9 Li systems, i.e. going vertically 
upwards in fig. [I] approaching the threshold, the binding decreases and the radius increases 
corresponding to a logarithmically diverging curve in fig. ^| (filled triangles) eventually ap- 
proaching that of the Borromean hypertriton example. A similar behaviour is seen for a 
hypothetical hypertriton with only s-waves included (filled squares) but with a slightly in- 
creased nucleon-A attraction. Then we move horizontally from the tango to the Borromean 
region. For a Borromean system of three particles (stars) where the three masses and scatter- 
ing lengths are equal or differ substantially, we again follow the same trajectory by changing 
one scattering length. 

By increasing the s-wave attraction in the neutron- 9 Li systems, i.e. vertically approaching 
the Efimov limit in fig. [j], larger binding and smaller radius result as expected. The resulting 
non-Borromean n Li ground state resembles more and more an ordinary nuclear state. How- 
ever, at some point an excited state appears, i.e. the first Efimov state (open triangles). At 
first it is very weakly bound and close to the Efimov line. As the attraction and the three-body 
binding energy increases the size at first decrease and then "turn around" and increase again 
as the binding energy of the two bound two-body subsystems approach and finally overtake 
the three-body binding energy |j| . As the two-body threshold is approached the third excited 
state (second Efimov state) should appear on the dashed line. Unfortunately the mass ratio 
(neutron to 9 Li) is relatively small and the next state is many orders of magnitude outside the 
scale of the plot, i.e. outside the reach of experimental as well as most numerical techniques. 

One remaining question is the approach of the two sets of points, related to hypertriton and 
11 Li, at very small binding energy. Both these sets arise by approach of (different) thresholds 
for Borromean binding in fig. [I]. When the majority of the radial wave function is located in 
the tail of the adiabatic potential in Eq.(^) then a av is a decisive length parameter. One could 
then erroneously be led to conclude that a av determines the size-binding relations and the 
threshold for Borromean binding in fig. [|. However, the absolute scale of the energy can not 
be determined from one parameter alone, the short distance behaviour is also indispensably 
necessary §,@,0. 

To understand this we imagine that we reduce the ranges of the potentials while increasing 
the strengths to maintain the scattering lengths. This zero-range limit results in the Thomas 
collapse and infinitely many bound states at small distances H]. To avoid this unphysical 
behaviour some kind of renormalization is needed. The simplest is to maintain the large 
distance behavior while only using the two-body potentials down to a distance below which 
the structure is uninteresting. The two-body results can then be expressed in units of such a 
rather arbitrary length parameter. The effect on a three-body system can then be anticipated 
mimicked by a similar renormalization by use of a hyperspherical length unit. However, this 
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is precisely the content of fig. where we used po as the scale parameter. 

The curves in fig. |^ coinciding at small binding are therefore an almost universal curve 
as indicated by the convergence in the figure. It would perhaps be rather fortuitous if this 
simple renormalization procedure in hyperradius results in a universal rescaled curve. The 
average scattering lengths for the four cases in the weak binding limit vary from 4 fm to 
20 fm whereas the ratio a av /po varies between 1.7 and 4.2. However, a av only determines 
the adiabatic potential for distances smaller than the scattering lengths and larger than the 
potential ranges. These conditions are not well fulfilled for the examples in fig. || Thus 
constant a av / p should not necessarily arise, since a av should be replaced by a complicated 
function of all three a^. This is not contradicting the numerical results in fig. |^ which is 
obtained without use of a av . Thus the emerging numerical curve is almost universal providing 
rather well defined scaling properties. 

Summary and conclusion. Three-body quantum halos were previously believed to 
appear along the K = curve This conclusion was reached by omitting constant terms 
compared to the leading order logarithmically diverging term. However, the present more 
refined analysis reveals the correct higher lying universal curve approached in the weak binding 
limit. In fact asymptotically these curves differ by a constant factor. In any case three-body 
quantum halos can appear above the K = curve, i.e. allowed in a larger window and being 
a factor of 2-3 larger than expected. One implication is that K = wavefunctions cannot 
be used in the weak binding limit due to the coherent large-distance contributions from more 
than one subsystem. For stronger binding a K = basis is simply incomplete. 

In fig. [I] ground state three-body halos appear just below the line separating the Borromean 
from the unbound region. On top of this the Efimov effect can also give "super-halos" in 
excited states. Ground state two-body halos can appear just below the line separating the 
tango from the unbound region. In fig. ^, three-body halos where the three-body dynamics 
dominate will appear on or below the (almost) universal curve. All systems falling above are 
influenced by an effective two-body threshold and correspond either to the Efimov line (caused 
by the threshold in two subsystems) or the tango-unbound line in fig. [l| We stress that all 
quantities used to place a system in fig. |^ in principle can be measured experimentally and 
that our results therefore can be used to classify the three-body nature of realistic systems. 
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